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We calculate the universal pressure constants of a stack of A'^ membranes between walls by strong- 
coupling theory. The results are in very good agreement with values from Monte-Carlo simulations. 
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Membranes formed by lipid bilayers are important biophysical systems occuring as boundaries of organells and 
vesicles. In equilibrium with a reservoir of molecules, tension vanishes and the shape is governed by extrinsic curvature 
energy Eq- If a stack of membranes is placed between two parallel walls, violent thermal out-of-plane fluctuations 
of the membranes exert a pressure p upon the walls. The pressure law was found by Helfrich Q| and reads for N 
membranes 
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where L = (N + l)a is the distance between the walls, and k the bending stiffness. The universal pressure constants 
aN are not calculable exactly. For a single membrane, ai was roughly estimated by theoretical [Q| and Monte-Carlo 
methods H-^. The most precise theoretical result was obtained by strong-coupling theory ^ yielding af^ — 0.0797149 
which lies well within the error bounds of the latest Monte-Carlo estimate af^ = 0.0798 ± 0.0003 [|. 

For more than one membrane between the walls, the strong-coupling calculation of Ref. ||^ must be modified in 
a nontrivial way. We must find a different potential which keeps the membranes apart and whose strong-coupling 
limit ensures non-interpenetration. For this we take advantage of the fact that membranes between walls have similar 
properties as a stack of nearly parallel strings fluctuating in a plane between line- like walls , in particular the same 
type of pressure law (1.1) with k substituted by the string tension a. The characteristic universal constants of the 
latter system were exactly calculated in Refs. |^,^ from an analogy to a gas of fermions in l-f 1 dimensions [|-[ll). We 
use these exact values to determine a potential which, when applied to the stack of membranes, yields a perturbation 
expansion for the pressure constants for an arbitrary number of membranes to be evaluated in the strong-coupling 
limit of complete repulsion. 

Our results are in excellent agreement with all available Monte-Carlo estimates for N — 1,3,5. By an 

extrapolation to A'^ — > 00 we determine the pressure constant aoo for infinitely many membranes. 

II. STACK OF STRINGS 

We begin by studying the exactly solvable statistical properties of a stack of N almost parallel strings in a plane, 
which are not allowed to cross each other and whose average spacing at low temperature is a. The system is enclosed 
between parallel line-like walls with a separation L as illustrated in Fig. ^ In the Monge parametrization, the vertical 
position of a point of the nth string is z„ = z„(a;). Since the vertical positions of the nth string are fluctuating around 
the low-temperature equilibrium position at na, it is useful to introduce the displacement fields ipnix) 



ipn{x) = Zn{x) - na. 
The thermodynamic partition function is given by the functional integral 
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where a is the string tension, T is the temperature, and fc^ the Boltzmann factor. We are interested in the free energy 
per unit length 
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with A — dx. Since the strings may not pass through each other, the fluctuations ipn{x) of the nth string are 
restricted to the interval 



Vn{x) e {ipn-l{x) + a, (pn+l{x) - o} . 



(2.4) 



A. Free Fermion Model 



The restriction (2^) makes it difficult to solve the functional integral {2/2) explicitly. It is, however, possible to find 
a solution using an alternative of the same systems as a 1 + 1-dimensional Fermi gas observed by de Gennes . Using 
this analogy, Gompper and KroU Q determined the 1/a^-contribution to A/^ relevant for the pressure law (1.1) as 



with the pressure constants 
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FIG. 1. Stack of A'' strings with equilibrium spacing a between two walls of distance L. The magnifier shows the local 
displacement field (pn{x) as the distance from the position Na. The walls are labeled by and A^+1 and treated as nonfluctuating 
strings with (po{x) = ipN+i{x) = 0. 

The analogy with fermions cannot be used to calculate the free energy of a stack of membranes, where only 
approximate methods are available. We shall use a strong-coupling theory as in Ref. ||^. As a preparation, we we 
apply this theory to the exactly solvable system of a stack of strings. 



B. Perturbative Approach 



The difficulty in solving the functional integral (2.2) arises from the restriction ( p.4| ) of the fluctuations by the 
neighbouring strings. To deal with this strong repulsion, we introduce into the action of the functional integral (2.2) 
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a smooth potential which keeps the strings apart in such a way that the integration interval for th e fluctuations can 



be extended to ipn{x) € {—00,00}. At the end we take a strong-coupling limit which ensures (2A)- In Ref. ||6|, such 
a method was used to evaluate the pressure constant for one membrane between walls. The smooth potential for the 
analogous case of one string is V{ip{x)) = {2a ^/i:Y tan^[7r(/?(a;)/2a] which describes the hard walls exactly for /i — > 0. 
This potential is symmetric and possesses a minimum at Lp{x) = 0. Thus its Taylor expansion around this minimum 
is a series in even powers of (p{x). 

In the case of N strings, the minima of the repulsion potential should lie at the equlibrium positions of the strings. 
The Taylor expansion of such a potential will also have terms with odd powers. Unlike the one-string system, where 
fluctuations are limited by fixed walls, the range of the displacements ipnix) of the nth string in an iV-string system 
depends on the positions z„_i(x) and Zn+i{x) of the neighbouring strings. Thus the potential will be taken as a sum 
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where WnZnix) denotes the pre-point lattice gradient Zn{x) — z„_i(a;). This potential includes the interaction of the 
first and last strings with the walls as non-fluctuating strings at zo = and z^+i = {N + l)a = L: 

ipo{x) = (fN+iix) = 0. (2.8) 

In the limit ^ — > 0, the potential V^(V„z„(a:)) should again yield an infinitely strong repulsion of two neighbouring 
strings for Zn{x) close to Zn-i{x). For z„(a;) > Zn-i{x), the limiting potential should be zero. 




FIG. 2. Potential V^(VnZn(a;)) of string-string interaction for finite /x and small /i as function of V„2:„ (a;) — z„{x) — z„-i{x) . 
The strings repel each other strongly for V„Zn(x) 0, while the potential has a minimuim at the equilibrium separation 
\/nZ„{x) — a, and we choose to normalize it to zero at that point. 

As a matter of choice, we let the potential between two strings V^(V„z„(a;)) be minimal and zero at the positions 

= na and z^l^ = a{n ~ 1): dF^(a)/d V„2„(x) = and F^z^^ - z^i) ^(a) = (see Fig. |). 
The Taylor expansion aromid the mimimun is, in terms of the variables (2.1), 



F^(V„(p„(x)) = ^[V„ip{x)f+f,^ J2 Cfc/[V„v'(x)]^+2. 
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The parameter fi governs the harmonic term. The introduction of the coupling constant g = l/a makes the coefficients 
Ck dimensionless, and the partition function (2.2) becomes 
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with the integral measure 
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The harmonic part of the partition function can be written as 
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with the functional matrix 
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Here V„(/3ri(x) = (pn+i{x) — (Pn(x) denotes the post-point lattice gradient in the 2;-direction, and V„V„ is the lattice 
version of the Laplace operator Q . 

Let us now impose the vanishing of the fluctuations of the walls ( |2.8| ) , corresponding to Dirichlet boundary condi- 
tions. For a finite number N of strings, the Kronecker symbol dnn' in Eq. (2.13) has the Fourier representation 
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with wave numbers Vm = TTm/{N + l)a. Thus the kernel [G^i„/ (a;, x')] ^ may be written in Fourier space as 
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Integrating over k in the spectral representation (2.15) leads immediately to the correlation function in configuration 
space 
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and to the harmonic partition function ( p.l2D 

Z;; = exp|-iTrln[G^]-i| = e-'^^«.f'/'=^^, 

the exponent giving the free energy per length: 
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which vanishes for /x = 0. 

The full partition function in Eq. (2.1C) is now calculated perturbatively. We introduce harmonic expectation 
values 
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in terms of which the correlation function is given by 

The perturbation expansion contains the two-point correlation function of V„iy9„(a;) which is given by 
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We now expand the second exponential in Eq. ( 2.10| ) in powers of the coupling constant g. Harmonic expectation 
values with odd powers of V„(p„(x) do not contribute, and the expansion reads: 
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In the sequel, we r estrict ourselves to the terms of second order in g — l/a, which contribute directly to the pressure 
law as in Eq. (2.5). The higher powers diverge for /i — > 0, and in Ref. p] it was shown how to calculate from them 
a finite strong-coupling limit. Here we shall ignore these terms for reasons to be explained shortly. Reexpressing the 
right-hand side of Eq. (2.23) as an exponential of a cumulant expansion, we obtain the free energy per length 
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We have used that the free energy ^ of the harmonic system ( 2.19 ) vanishes in the limit fi 0. The first cumulants 
are the expectations 

(Oi(V^„,(a;i)));_^ - (Oi(V(p„,(xi))); 

{Oi{W^m{xi))02(W^nAx2)))l^^ = (Oi(Vv7„,(xi))02(V(^„.(a:2)));,-(Oi(V(^„,(a;i)))^^(02(V(^„.(a;2)))^, 

: , (2.25) 

defined for any polynomial function Oi(V(/J„; (x^)) of 'S/ipniixi)- Following Wick's rule, we expand the expectations 
into products of two-point correlation functions (2.22). The different terms are displayed with the help of Feynman 
diagrams, in which lines and vertices represent the correlation functions and interactions: 

xi,ni X2,n2 > (V„,(p„,(xi)V„2(^„,(x2))^^, 
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In what follows, we assume that the potential parameters Ck with k > 3 are chosen in such a way that they make all 
terms of order g^ and higher vanish. Dividing the free energy (2.24) by N, we obtain the following expression for the 
free energy per length and string which can be compared with Eq. (2.5): 
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The calculation of the Feynman diagrams is straightforward using Eq. (2.17). The evaluation is only complicated by 
the Dirichlet boundary conditions, which destroy momentum conservation. This makes the numeric calculation quite 
time-consuming for an increasing number N of strings. As an explicit example consider the sunset diagram which 
requires the evaluation of the multiple sum 
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with the abbreviation 
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It is useful to factor out the physical dimension of the diagram. Any Feynman integral with I lines and v vertices 
can be expressed in terms of a reduced dimcnsionless Feynman integral W^^'^ as 
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This brings Eq. ( ^.28| ) to the form 
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where the diagrams indicate the reduced Feynman integrals. Their values are listed in Table ^ for different string 
numbers N. Note that the l/a^-contributions to the free energy per length and string in Eq. (2.28) is independent 
of /i since the /i-prefactors are cancelled by the /i-dependence of the diagrams. Thus the limit fi ~f becomes trivial 

for these contributions. 

With the knowledge of the exact values of the c onstants a% from Eq. ( |2.6| ), we are now in a position to dete rmin e 
the potential par amete rs ci and C2 from Eq. ( ^.33 ) to obtain the exact result from the two- loop expansion Eq. ( 2.33 ). 
Comparing Eqs. (2.33) and (p^) for = 1 and A^ = 2, we obtain 



Cl 



77 

3' 



C2 = 



6 



(2.34) 



Note that (2.33) consists of more equations than necessary t o co mpute ci and C2. It turns out, however, that all of 



them give the same ci and C2, such that the same potential (2.9) can be used for any A^. This is the essential basis 
for applying this procedure to a stack of membranes. 

We now justify the neglect of the higher g powers which would in principle give a further contribution to the 
pressure constant in the strong-coupling limit. We simply observe that it is possible to choose the higher expansion 
coefficients to make all higher contributions vanish . 



III. STACK OF MEMBRANES 



Having determined the parameters Ci and C2 of the Taylor expansion (|]^) of the smooth potential applicable for 
any number of strings, we shall now use the same potential for a perturbative expansion in a stack of A^ membranes 
displayed in Fig. |3[ 
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FIG. 3. Stack of selfavoiding fluid membranes fluctuating in ^-direction between two walls. As for the previous stack of 
strings, the walls are treated as nonfluctuating membranes. 



The equilibrium spacing at low temperature between the membranes is again a. Denoting the vectors in the plane 
by X = ix,y), the vertical displacements of the membranes from the positions na are (p„(x), with Dirichlet boundary 
conditions at zq and zj^^i 



For membranes without tension, the energy has the harmonic approximation 
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where k is the bending stiffness and = + dy is the Laplacian in the plane parallel to the walls. By analogy with 
the previous section, the kernel of the harmonic stack reads now 
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where we have used a mass parameter /x* instead of /i^, for dimensional reasons. The partition function for the stack 
of membranes is then written up to order = by 
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with the same parameters ci and C2 of the Taylor expansion (2.£) as in the string system, determined in Eq. ( ^.34 ). 
We neglect terms of order which certainly contribute in the strong-coupling limit, and which vanish only for the 
strings, where the partition function ( p. 23 ) with the choice (2.34) for the parameters ci, C2 is exact in second order. An 
evaluation of the neglected terms by variational perturbation theory is expected to give only a negligible contribution 
to our final result. 



Inverting the kernel (3.3) yields the correlation function 
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The explicit calculation of the Fourier integral leads to a difference of modified Bessel functions -^^0(2;) as in Ref. 
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For Xi = X2 = X and ni = 712 = n, this reduces to 
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leading to the partition function of the harmonic system 
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where A — ^ cPx is the area of the projected plane of the membranes. The free energy per area fjy^^ — — {kBT/A)ln 
vanishes again for ^ = 0. 

As for the calculation of the free energy of the stack of strings, we introduce harmonic expectation values 
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which appear in the perturbation expansion of (3.4), the cumulants yielding a perturbative expansion for the free 
energy per area Jn = —{kBT/A)lnZ. The lines and vertices in the Feynman diagrams stand now for 
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and the two-loop approximation to the free energy per area and membrane in order l/a^ reads 
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Going over to reduced Feynman integrals as in Eq. (2.31) 



where v is the number of vertices and I the number of lines of the diagram, we obtain 
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The pressure exerted by the membranes upon the walls is obtained by differentiating the free energy /jv = NAf^ 
with respect to the distance of the walls L = a{N + 1): 
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The first and the last Feynman integrals in Eq. ( 3.15| ) are the simplest: 
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where we have used the abbreviation /i™^^ defined in Eq. (2.30). The evaluation of the second diagram in Eq. ( 3.15| ) 



is much more involved. The Fourier integrals can be done exactly, except for one, which must be treated numerically. 
This calculation is deferred to Appendix The values of the three diagrams are listed in Table for various numbers 



of membranes. With these numbers, the evaluation of the pressure constants yields the results given in Table II] 
Except for = 1 and N oo, no analytical values were found in the literature. We also compare with pressure 
constants obtained by Monte-Carlo simulations and find a good agreement [^,1). The values of the Monte-Carlo 
simulations for A^ = 3, 5, 7 from Ref. |^ show an independence of the number A^ of membranes. This arises by the 
simulation technique, where the free energy of the central membrane was determined. In contrast to that, we have 
calculated the pressure constant from the free energy of the complete system averaged over all membranes. Thus 
these Monte-Carlo values cannot be directly compared with ours. 



Table III contains also a value for an infinite number A^ ^ oo of membranes in the stack. This pressure 
constant is obtained by the following extrapolation procedure. We assume that the pressure constants determined for 
A^ = 12, 13, 14, 15 are of higher accuracy than those for lower numbers of membranes. This assumption is justified by 
comparing our values for A^ = 1, 3, 5 with the latest Monte-Carlo results For A^ = 1, the deviation is about 3.4%. 
Considering A^ = 3, the deviation reduces to 1.8% and further to 1.1% for 5 membranes. Since the pressure constants 
is approximated increasingly fast with increasing number A" of membranes, we make the following exponential ansatz 
for determining the approach to infinite A: 

ajv = aoo[l-77exp(-e^^")]. (3.19) 

The unknown four parameters in this equation are then determined by solving the system of equations with the 



pressure constants au, ai3, ai4, and ai5 listed in Table III. We obtain r/ sa 1.1712, ^ w 1.6417, e ~ 0.3154, and thus 



the limiting pressure constant for an infinite stack of membranes 

aoo ~ 0.1041. (3.20) 

This value is in very good agreement with the Monte-Carlo result |^ (see the last row of Table |l|). It differs by a 
factor close to 9/4 from the initial result by Helfrich |p],|l5|. 



IV. SUMMARY 



We have calculated the pressure constants for a stack of different numbers of membranes between two walls in 
excellent agreement with results from Monte-Carlo simulations. The requirement that the membranes can not pen- 
etrate each other was accounted for by introducing a repulsive potential and going to the strong-coupling limit of 
hard repulsion. We have used the similarity of the membrane system to a stack of strings enclosed by line-like walls, 
which is exactly solvable, to determine the potential parameters in such a way that the two-loop result is exact. 
This minimizes the neglected terms in the variational perturbation expansion, when applying the same potential to 
membranes. 
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APPENDIX A: EVALUATION OF THE SUNSET DIAGRAM 



The second diagram in Eq. ( 3.15|) requires some simplification before the numerical calculation. We write the 
reduced Feynman integral as: 
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All integrals are easily calculated, except for one. If we introduce abbreviations 
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we are left with integrals of the type 
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where 71^2 = iMi_2 are real numbers. Employing Feynman's parametrization, these integrals become 
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b + 2ax 



arctanz(x) 



6 = fc^ + z(7i - 72), c = i72. 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 



(AlO) 



Since b is a complex number. Re arctan z is discontinuous, if Re z changes sign and |Im z| > 1. Thus the right-hand side 
of Eq. (A9) is discontinuous at a certain point xq within the interval [0, 1]. As will be seen s ubse quently, J(k, Mi, M|) 
from Eq. (A5) must be real and thus all imaginary contributions in the decomposed form (A6) cancel each other. 

We determine the point of discontinuity xo to obtain the solution of the integral (A8) by investigating the zero of 
the real part of z{x). Decomposing z{xo) into real and imaginary part, we obtain 



Re2(x) = \A\-^/^ 
lmz{x) = |Api/2 



fc2(l - 2x) cos 
(71 - 72) cos 



1 ReA 
— arctan 

2 ImA 

1 ReA 
— arctan — 

2 ImA 



, . , 1 ReA 

+ (71 - 72) sm ( - arctan j—^ 

.1 R-eA 

— K (1 — 22;) sm - arctan - 

^ ^ * 2 ImA 



(All) 
(A12) 



where 



Re A = (71 - 72)2 - fc^ Im A = -2^^(71 + 72). 



(A13) 
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Thus, the zero of Kez{x) is found at 



Xo 



1 H "2 — tan 



Cretan ^^^(^1+^^) 



(AM) 



From the bounds of integration in Eq. ( |A8| ) follows that we must include the disco ntin uities of (A£) for xq e [0, 1]. 
This occurs if A: < I71 — 72I and sign 71 7^ sign 72. Thus the solution of the integral (A8) reads 



/(k,7i,72) = 



S'(k,7i,72,a::) 



lim 



x=0 

5'(k,7i,72,a;) 



sign 71 = sign 72 V (sign 71 ^ sign 72 A fc > ^/l-fi - 72I), 



x—XQ — e 
x=0 



S'(k,7i,72,a;) 



x—XQ+e 



sign 71 sign 72 A fc < \/|7i - 72I, 



where S'(k, 71, 72, x) is the explicit right-hand side of Eq. (A9): 

1 



S'(k,7i,72,x) 



27r V(7i - 72)2 - fc^ - 2ifc2(7i + 72) 
The function /(k, 71,72) possesses the properties 



: arctan ■ 



fc2(l-2a:)+z(7i-72) 

ye^r^^;F~fc^"^2iF(^rT^ ■ 



/(k, 71,-72) + /(k, -71,72) = 2Re/(k,±7i,=F72), 
/(k, 71,72) +/(k, -71,-72) = 2Re/(k,±7i,±72). 



(A15) 



(A16) 



(A17) 
(A18) 



Inserting E q (Al^ ) into Eq. ( |A6|) and performing the remaining integral in Eq. (A4) together with the sums in 
expression (Al) for the sunset diagram can be calculated numerically. The values are listed for = 1, . . . , 15 in the 
third column of Table ||. 
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TABLE I. Reduced numeric values VV'^ of the two-loop diagrams for the free energy for a stack of A'^ strings. 




11/2 

2 1.288675 0.398717 0.089316 

3 2.100656 0.832299 0.146447 

4 2.915827 1.270787 0.184463 

5 3.730993 1.709326 0.211325 

6 4.545586 2.147034 0.231245 

7 5.359574 2.583849 0.246583 



TABLE IL 


Numeric values W of the reduced two-loop Feynman integrals contributing 


to the pressure constants of a stack 


of N membranes in Eq. ([3.15[) . 








N 




C»' 






1 




1/32 








2 




0.080542 


0.022446 


0.005582 


3 




0.131291 


0.046992 


0.009153 


4 




0.182239 


0.071866 


0.011529 


5 




0.233187 


0.096762 


0.013208 


6 




0.284099 


0.121619 


0.014453 


7 




0.334973 


0.146428 


0.015411 


8 




0.385815 


0.171195 


0.016172 


9 




0.436630 


0.195925 


0.016789 


10 




0.487422 


0.220624 


0.017300 


11 




0.538197 


0.245300 


0.017730 


12 




0.588958 


0.269954 


0.018097 


13 




0.639706 


0.294592 


0.018414 


14 




0.690444 


0.319215 


0.018690 


15 




0.741174 


0.343827 


0.018933 



TABLE III. Pressure constants ajv for different numbers N of membranes in the stack, calculated from Eq. (3.15), with 
the numerical values of the two-loop diagrams given in Table We compare with results from Monte-Carlo simulations and 
earlier analytic results. 

Monte Carlo results earlier analytic values 

0.060 ig, 0.078 ±0.001 [|, 0.0798 ± 0.0003 §, 0.080 § 7r7l28 0.079715 |] 

0.093 ± 0.004 [|, 0.1002 ± 0.0006 § 

0.0966 1], 0.1022 ±0.0006 § 

0.1009 ±0.0007 [i 



N 


ajv 


1 


7r7l28 ^ 0.07711 


2 


0.08669 


3 


0.09134 


4 


0.09408 


5 


0.09590 


6 


0.09719 


7 


0.09815 


8 


0.09890 


9 


0.09950 


10 


0.09999 


11 


0.10039 


12 


0.10074 


13 


0.10103 


14 


0.10129 


15 


0.10151 


oo 


0.10409 



0.074 ||l|, 0.101 ± 0.002 §, 0.106 [| 37rVl28 fa 0.23 § 
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